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We use heavy quark effective theory (HQET) and perturbative QCD to 

> : 

' study the heavy meson - light vector meson transitions involved in D and B 

' decays. HQET is used to relate the measured D K*iv vector and axial- 

^ ! 

en ' vector form factors at four- momentum transfer = Qto the B K*^ tensor 

• and axial-tensor form factors at = 16.5 GeV^. Perturbative QCD is then 

Q^' used to find matching conditions for the i?-meson form factors at = 0. A 

i-S^ ' five parameter "vector dominance" type fit of the two HQET form factors, 



X 



consisting of single pole, double pole, and subtraction terms, is used to match 



■ the data at q^ = 16.5 GeV^ to QCD at q^ = 0. The values at = are 

compared with recent data on the exclusive rate for B K*j decay to extract 

a value for the Cabibbo-Kobayashi-Maskawa matrix element \Vts\ = 0.035, 

with 28% experimental uncertainty, and 32% theoretical uncertainty from 

higher order QCD effects and violations of heavy quark symmetry. 
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I. INTRODUCTION 



Rare B decays are a sensitive probe of flavor clianging processes. In particular, the 
quark level process b ^ s'~f occurs mainly via a "penguin" type diagram that is sensitive 
to the Cabibbo-Kobayashi-Maskawa matrix element \Vts\- A proper analysis of the penguin 
requires a treatment of leading logarithmic strong interaction corrections [|T| . At low energies, 
assuming only standard model physics, it corresponds to the effective Hamiltonian 

H^ = 7]mhSa^"'{l + -f5)bFf,^ + ... , (la) 

where r] = ^^^^^ Vtb V* F2 {mil m^) . (lb) 

The parameter F2 includes the leading logarithmic corrections and is mildly sensitive to the 
top quark mass. It ranges from 0.59 for = 120 GeV to 0.68 for mt = 210 GeV. We will 
use the value F2 = 0.63, which corresponds to rrit = 150 GeV. 

Recently, experimental evidence for this penguin has been found by the CLEO collabo- 
ration via the exclusive decay B — > K*j 0]: 

B.R.(5 ^ /s:*7) = (4.5± 1.7) X 10^^ . (2) 

This is related to the decay rate via F = B.R./r^, where we use the B meson lifetime 
tb = (1.29 ± 0.05) ps The theoretical expression for the exclusive rate of is given by 

where A^^\{}i) and i?(^^(0) are tensor form factors @] determined by the matrix element 
{K*\s + 75) h\B) at = and defined explicitly in section two. In this paper, we will 
study the consistency of the standard model with this data, by extracting from it a value 
for \Vts\- Throughout this paper we will take Vtb = 1 in the context of the standard model 
with three generations 0, although it should be clear from eqn. (lb) that we are really 
calculating the combination |VJs^6|. 
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Unitary of the Cabibbo-Kobayashi-Maskawa mixing with three generations of quarks 
places a severe constraint the mixing angle \Vts\- Following the conventions of the particle 
data group 0, \Vts\ ~ = 0.041±0.007. An additional quark generation or the pres- 

ence of vectorlike down type quarks (which appear in SO{10) unification) would invalidate 
this result. The penguin itself is quite sensitive to non-standard model physics, particularly 
to the exchange of charged Higgs bosons which can easily suppress or enhance the 

parameter F2 by a factor of two. In principle, this can be used to place constraints on 
supersymmetry breaking parameters 0. 

Quark model calculations of the hadron matrix elements and the subsequent branching 
ratio vary by an order of magnitude and are therefore unreliable as a test of the standard 
model result for \Vts\- Instead, we shall apply heavy quark effective theory (HQET). This 
is a model independent framework that relates processes where a heavy quark of mass 
rriQ ^ Aqcd exchanges momenta smaller than mg with the light degrees of freedom inside 
the hadron. HQET makes manifest the symmetries of the QCD Lagrangian that occur in the 
infinite heavy quark limit: a flavor symmetry which relates the mass splittings and decay 
amplitudes of hadrons with different heavy quark content, and a spin symmetry which 
simplifies the mass spectrum and relates the decay amplitudes of hadrons with the same 
heavy quark content, but with different heavy quark spin. The spin symmetry is analogous 
to the proton spin symmetry in the hydrogen atom spectrum, which is broken only weakly 
by the hyperfine splittings. (For a review of HQET, see and references therein.) For 



example, all of the form factors relevant in the semileptonic decay of a i? to a -D or D* 
(charmed) meson can be given in terms of a single universal function, to leading order in 
an expansion in Kqcd/^q-, where mg corresponds to h and c quarks masses. This so called 
Isgur-Wise function depends only on the relative four velocities of the initial and final heavy 
quarks, and is absolutely normalized at zero recoil. 

HQET also has implications in heavy hadron to light hadron (heavy - light) transitions. 
The process B -rIv ,i.e., the decay of a heavy pseudoscalar meson {s'^ = 0") to a light 
pseudoscalar meson, has recently been studied in this context [jlll. For the case at hand. 



{B — > K*'y), we will need to consider the decay of a heavy pseudo-scalar meson to a light 
vector meson [s'^ = 1^). Therefore, in Section 2 we will use the tensor formalism |I2[ to 



establish the most general form of any transition amplitude between a heavy pseudoscalar 
meson and a light vector meson consistent with HQET. Again, the strange quark is not 
treated as heavy compared to the QCD scale Aqcd- 

To extract \Vts\ from the data, we need to determine the combination of tensor form fac- 
tors A('')(0) + fi('')(0) which occur in the B K*'y decay rate, eqn. (|^). Here we apply heavy 
quark symmetry to relate this combination to the form factors associated with D K*ii' 
dataQ, 

{K*ip', e)\sr{l - 75) c\Dip)) = elp'^Pf^Viq') - (m^ + mK*)e*^A,{q') 

+ e* ■ p{p + p'Y + — M^e* ■ p{p - p'Y , (4) 

niD + rriK* rrif) + rriK* 

where = {p ~ p'Y- The form factor A3 leads to unmeasurable contributions proportional 
to the electron mass. For Ai, A2, and V, we will use the averages over data from experiments 
E653 {D K*ev) and E691 ^ {D ^ K*fiu) at Fermilab given in refs. 



1/(0) = 0.95 ± 0.20 , Ai(0) = 0.48 ±0.05 , ^2(0) = 0.27 ± 0.11 . (5) 

These are the values of the form factors given at zero invariant momentum transfer, = 0, 
for the D —>■ K*lv system, extrapolated with a simple pole ansatz for the form factors. (The 
results are relatively insensitive to the form of the ansatz, because the data is taken in the 
narrow region = to ~ 1 GeV^.) 

Heavy quark symmetry will relate the two decay systems in in the following way: since to 
leading order the heavy quark form factors calculated from the effective theory are indepen- 
dent of the heavy quark mass, they are functions only of the K* mass and the dimensionless 
parameter 



We use the convention eoi23 = +1- 
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w = , (6) 

rriK- 

where v is the velocity of the heavy meson and p' is the momentum of the final state K* 
meson. On the other hand form factors are given as a function of 

= tuq + - 2mQmK*w , (7) 

where itlq is the mass of the heavy meson. This means that = for the B K*'y decay 
corresponds to w = wb ~ 3.04, while for the D K*ii> decay, = corresponds to 
w = wd ~ 1.29. We will show in section 2 that the heavy quark flavor and spin symmetries 
relate the D K*lv data given by (|^) to the tensor form factors of the B K*'y decay 
at w = Wf), which corresponds via eqn. (0) to ~ 16.5 GeV^ for the B decay system. We 
must then evolve the B decay form factors down to = 0, or equivalently, from m7£) up to 
Wb- 

To determine the interpolating functions that can be used to evolve the form factors from 
Wd = 1-29 to w = 3.04, which for B K*'y decay is the kinematical point corresponding to 
the emission of an on shell photon, we will use Perturbative QCD. This is the appropriate 
description of strong processes involving the exchange of hard gluons, and corresponds to 
heavy - light transitions with a large value of the kinematical variable w, that is, far away 
from the "zero recoil" point w = 1. Therefore, in Section 3 we develop the Brodsky-Lepage 



formalism |]T^ of perturbative QCD to calculate the HQET heavy - light form factors for 
"large" w, to leading order in the heavy quark mass expansion, and to leading order in as- 
We argue that the method is reliable at w = wb, (corresponding to ~ 0.20), and use this 
perturbative calculation to place "matching constraints" on the interpolating form factors. 
Furthermore, we use perturbation theory to determine the leading violations of HQET - the 
order w ■ mx* / fng corrections. 

We will use the D K*[v data to constrain the interpolating form factors at Wd-, and 
perturbative QCD to constrain them at wb- As discussed in section 4, this will enable us 
to make a five parameter fit to the two relevant heavy-light form factors required for the 



determination of \Vts\. The parametrization is in terms of single pole, double pole, and a 
single subtraction (constant piece), consistent with "vector dominance" ideas. 

This is a better way of fitting to perturbative QCD than simply assuming a pole form for 
the vector and axial-vector form factors because we require the parametrization to match 
the heavy - light form factors in the perturbative regime. This is not the same as matching 
the perturbative result in the limit w oo, with rriQ fixed; in this limit, perturbative 



QCD counting rules [|T8l do indeed indicate single pole dominance, but this is due to the 
w ■ mK*/TnQ corrections from the point of view of the HQET. The HQET limit corresponds 
to taking rriQ —>■ oo, with w fixed, and this is the correct way of performing perturbation 
theory in our case, since the evolution of the form factors occurs for w ■ mK*/'mB < 1. In 
fact, as we shall show in section 2, assuming a simple pole form for these form factors, with 
no subtactions, is inconsistent with heavy quark symmetry to leading order in heavy quark 
mass. 

As a self consistency check, we use the absolute normalization of the B decay form factors, 
obtained by HQET matching to the D decay data, to estimate of the B-decay constant 
to leading order in QCD perturbation theory. This check, and further error estimates, are 
contained in section 5. 

Section 6 is devoted to a discussion on how to systematically improve these results, and 
our conclusions. 



II. HEAVY LIGHT FORM FACTORS 

We shall first apply the tensor method [|1^ to the case of a heavy pseudoscalar meson 
Mqiy) decaying into the light vector meson K*{e,p'). This method is usually applied to 
the derivation of heavy to heavy meson form factors P,p!0| , and recently has been used for 
the case of the heavy meson to light meson transition B niu |[Tl]|. We will determine the 
most general form of heavy psuedoscalar to light vector form factors consistent with Lorentz 
invariance and heavy quark spin symmetry. 



6 



It is useful to consider the matrix element of the general operator sTQ, where T is an 
arbitrary product of gamma matrices, and Q is a heavy quark Dirac field. To leading order 
in the HQET, Q is replaced by its projection onto the quark (versus antiquark) field via the 
replacement Q h^^\ where h^\x) = exp {iniQV ■ x)^^Q{x), and v is the heavy meson 
velocity obeying = 1. The relevant matrix element takes the form 

(ir*(p^e)|s^/l(Q)|MQ(^;)) = y^7^F^Tr|e(^;,p^e*)^ii^75| , (8) 

where p'^ = mK*v' ^ is the K* momentum, and e is the K* polarization vector which satisfies 
e* ■ p' = 0. The in Eq. (H) explicitly projects the matrix element onto the heavy quark 
components of the heavy quark spinor, since in the infinite quark mass limit heavy antiquarks 
are not produced, and the 75 describes the pseudoscalar nature of the heavy meson. A simple 
way to understand eqn. (Q) is by parametrizing the initial heavy meson as j/i^'^Hs^), where 
C. describes all of the light degrees of freedom in the heavy meson, and using the heavy 
quark propagator (h^^'^h^^) = (l + ^)/2 . The matrix function G encapsulates all of our 
ignorance of the dynamics of the light degrees of freedom of the heavy meson and the K* 
meson. It is explicitly independent of the heavy quark mass to leading order in HQET. The 
y/mq in eqn. (|^) comes from the usual normalization of the heavy meson with respect to its 
energy. 

More specifically, the scalar matrix B must be proportional to the polarization of the K* 
meson: G = Aii^* + ■ e*, where the Aii are matrix functions of rriK* and w defined in 
eqn. (|^). Furthermore, using the fact that = the only possible matrix structure of 
each Aii within the trace is + f'9\ where the 6 and 0' are real functions of v-p' . Therefore, 
one can parametrize B in terms of only four linearly independent form factors, 6i{mx*,w), 

e = (01 + ij'e2)f + (e-i + ij'e^)v ■ e* . (9) 

We will refer to Eqs. (§|J^) as the heavy - light matrix elements for the decay of a heavy 
pseudoscalar into a light vector. This is the relativistic generalization (where the heavy 
quark is not necessarily in its rest frame) of the heavy - light matrix elements obtained in 



Since the b and c quarks are not infinitely heavy, the relation between operators in 
QCD and in the heavy quark effective theory is non-trivial. In the leading logarithmic 
approximation, the relation between the QCD currents and heavy quark effective theory 
currents discussed above is sTQ = CQ{fi)sTh[^\ where the coefficient functions are ||2^ 
Cq^ii) = [a{mQ) / a{jji)]~^^^'^'^~'^^\ N denotes the number of flavors below the uiq scale, 
and jjL is the infrared subtraction point. It is useful to define 6'|*^'* = CqOi. The 6*1*^^ are 
subtraction point independent and satisfy 



n{b) _ ^ q{c) ^ 



a{mb) 



airrir 



-6/25 



(10) 



We use as(mfe)/a^(mc) = [In {ni^/KQCD}/ In {m^/ Xqcd)], which for Kqcd = 300 MeV gives 
Ccb = 1.114. 

Eqs. (|^, place strong constraints between the form factors corresponding to different 
matrix elements. For F = a^'^, the matrix element is given by 

{K*ip',e)\rsa'^'^Q\MQip)) = e^^"^ (A(«)e>^ + 5(«)e>'^ + C^^^e* ■ . (11) 

Using the identity cr^'^'y^ = —^i€^'^°'^aa/3 one can immediately write the form factors for the 



F = matrix element. 



{K%p',e)\sa^^j,Q\MQ{p)) = z 



A(e*>^ - e* V) + B{e*^p"' - e* >"") 
+ Ce* -p {p^'p'" - pV^) 



(12) 



Evaluating these matrix elements using Eqs. (§J^) gives relations between the "tensor" form 
factors A, 5, C, and the heavy - light form factors 6'j, 



29 



(Q) 



3/2 1/2 • 

The vector and axial-vector matrix elements can be parametrized as 



(13) 



{K*{p',e)\srQ\MQ{p)) = zD^^^e^'^'^Pp^ey 



(14a) 



(ir* (p', e) I s 7^75 Q I Mq(p)) = E^^h*^' + F^^h* ■ pp^ + G^^h* ■ pp'^ . (14b) 
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Evaluating these matrix elements using the heavy - light parameterization yields 

DiQ)^^^^^ EiQ) ^ 2._frfu77m?]ef^ - - OP . (15a) 

dQ) n(niQ) , n(.Q)\ 

= -2V^7I^^ , = (15b) 

Since the seven form factors for the matrix elements are given in terms of only four heavy 
- light form factors, there are three relations relations between them for each heavy meson 
Mq: 

^^.EHp-p')D^ C = (16) 

These relations explicitly relate the tensor form factors to the vector and axial-vector form 
factors for a given heavy meson Mq transition to K*. They represent the explicit realization 
of the heavy quark spin symmetry for the heavy pseudoscalar meson decay to a light vector. 
Equivalent relations, determined via analysis in the heavy quark rest frame, are found in 
refs. PJTm. This confirms our analysis which led to eqns. (|^J^). 

The relationships between the heavy - light form factors at the kinematical point wo = 
(m|) +111]^,) /2m £,771 K* ~ 1-29 and the experimentally measured D — > K*lv form factors are 



771 £) + 771k* 
2^7nD77lK* 
y/77lK*r7lD 

7nK* + rriD 



MO) 
no), 



ml-* + 



{7171) + TTLk* 



771 K* ]] 77iK* mo + rriK* 



(17a) 
(17b) 
(17c) 
(17d) 



rrix* 771 D + 771k* 

From the data, the values of ^1, 62 and the combination 9^ — (m£)/mx*)6'4 at wd can be 
extracted. 



e^{^{wD) = -(0.06 ± 0.13) , ^(wd) = -(0.44 ± 0.09) 



(18a) 
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et\wD) - ^e^l^{wD) = -(0.66 ± 0.22) , (18b) 

where we have quoted experimental uncertainty. 

Treatment of the strange quark as heavy in the context of HQET has been considered in 
the hterature Eqns. (|1^) indicate that Aqcd/^s corrections are large. In an expansion 
in terms of the strange quark mass, the trace formula eqn. (|]) for the heavy - light form 
factor is modified by addition of ^*(1 + '^')/2, to project onto the "heavy" strange quark 
(verses antiquark). This exercise gives in the standard result for the Isgur-Wise function 
^. Using the heavy - light parametrization, one finds that ^ = 61 — 62, and further that 
61 + 62, 63, and 64 each vanish to leading order in the l/rris expansion. Clearly, the vanishing 
relations are strongly violated by the data given above. 

Heavy quark flavor symmetry relates the heavy - light form factors for D decay to those 
for B decay, as given in equation (|l^). The data determines 61 and 62, which in turn give 



the tensor form factors required for the extraction of \Vts\ via eqn. (13). We are then left 
to determine the running of the form factors from wd to wb- This is a dynamical question, 
and the kinematical heavy quark symmetry cannot give us the answer. 

There are two "straightforward" avenues of approach to this problem, QCD sum rules 



22 1, and perturbative QCD |0. In the next section, we will discuss the later approach. We 



note that in the literature, a third approach is discussed, which is to make an outright guess 
as to the dependence of the form factors, following vector dominance ideas. To see the 
danger in this, consider assuming the standard single pole vector dominance ansatz for the 
vector and axial vector form factors V{q) and Ai{q) for the B meson decay. This means that 
V and Ai are of the form m\/{q^ — m"^), where rriT a vector or axial vector threshold mass 
above the B meson mass. Via eqns. (^J^), in terms of the more relevant HQET variable w, 
single pole dominance is of the form 1/w, plus order Aqcd/^q corrections. However, now 
that we have correctly parametrized the heavy - light form factors, it becomes clear that 
this form is incompatible with leading order HQET. Since V oc 62, and Ai oc (^^i —W62), Ai 
must have a "large" subtraction if the vector form factor obeys single pole dominance. 
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III. B K*-f FORM FACTORS AT LARGE W FROM PERTURBATIVE QCD 

In the previous section we have used HQET to relate the form factors in D — K*ei' 
with those in B ^ K*'y at wd = 1.29. However the physical value of w for B K*'~f that 
corresponds to an on shell photon is wb = 3.05. This value for the B decay is in the regime 
of large w far from zero recoil. How does one estimate the behavior of the form factors in this 
regime? First consider the treatment of heavy meson to heavy meson transitions by heavy 
quark symmetry. All form factors are determined from the Isgur-Wise function C.i'w) and 
heavy quark symmetry fixes the value C,{w) = 1. For large enough values of w the Brodsky- 
Lepage method of determining form factors from perturbative QCD can be invoked. To 
leading order, the Brodsky- Lepage formalism essentially includes the exchange of a single 
hard gluon between the initial or final heavy quark and its spectator quark. Assumptions 
about the meson wave functions, motivated by QCD sum rules and experiment are made 



to fix the soft behavior of the amplitude. This method determines p3[ the behavior of the 
Isgur-Wise function at large w. For itlq/ Kqcd > w ^ 1, ^(w) ~ -^^^r^- As discussed in 
the previous section, this corresponds to a dipole in the sense of a vector dominance model. 
For the true asymptotic regime oi w > niQ/ Kqcd^ Ciw) ~ ^^^T^' "^^lere = —q^. This is 
in accordance with QCD dimensional- counting rules |]18| . 

We shall essentially follow the same procedure in discussing heavy meson to light meson 
transitions. In this case of course we have 61,62,03 and 64 instead of the single Isgur-Wise 
function. Although we don't have the convenient normalization at w = 1, we can use instead 
the D —>■ K*ev data to determine the thetas at w = 1.29. The fact that the kinetic point 
corresponding to B ^ K*'~f at w = 3.05 is far from zero recoil is a virtue; it is in just the 
right regime for both perturbative QCD and heavy quark symmetry to be valid. That is, if 
w is close to 1, perturbative QCD is not valid, while if w > mg/ Aqcd, then w ■ Aqcd/^q 
corrections dominate the HQET. 

We now briefly describe standard aspects of the Brodsky-Lepage formalism as it applies 
to the case at hand. The soft nonperturbative part of the physics is encapsulated in quark 
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FIG. 1. Two choices for the K* distribution amphtude. The sohd curve is (pi{x), and the 
dashed curve is (pczzi^^) j ^-s defined in the text. Both are normalized to unit area. 

distribution amplitudes (f){x,P,Q), for each meson, which denote the fraction < x < 1 of 
momentum P carried by the valence quark of the meson. To leading order in the calculation, 
the antiquark carries momentum (1 — x); it can be argued that gluon and non- valence quark 
corrections are small for large momentum exchange via the QCD dimensional-counting rules. 
The parameter Q denotes the subtraction point at which the distribution amplitudes are 
evaluated. Given a distribution amplitude (p at Qo, (piQ), for Q > Qq can be determined 
T%. The dependence is logarithmic in Q, and (p asymptotically approaches oc x(l — x). 



We shall initially consider two distribution amplitudes for the K* meson, as shown in fig. 1, 



^1 = Qx{l-x) . (19a) 



42 

-I2x^{l-xf + —x{l-x) . (19b) 



The second of these is the Chernyak-Zhitnitsky-Zhitnitsky distribution amplitude for 
the K* at = 1.5 GeV^ The integrals over X for the functions are both normalized. 
However, the correct normalization is given by = /x*0/2v^, where fx* is the meson 
decay constant (With this normalization, the pion decay constant /^r = 131 MeV 0). 
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For the momentum fraction of the heavy b quark inside the initial B, we use the peaking 
approximation. The distribution amphtude for heavy quarks must have a large maximum 



near the end point x = 1 — A/rriQ, and width A/mg, where mg = + A. We use |T0 
A = 0.5 GeV as the typical amount of mass in the heavy meson due to the light degrees of 
freedom and make the approximation 

M-) = ^^mi-^-^h (20) 

The error induced by this approximation is of order A/tuq, which we are systematically 
neglecting in this paper. 

The initial and final states are given by convoluting the distribution amplitudes with 
on-shell spinors of the quark and antiquark, 

^ = 0E/AA'WA(xP)t;A'((l-^)^) , (21) 

where A and A' denote sums over helicities. In the context of the peaking approximation, 
{If — mQ)u\{xP) = 0, and {If + mQ)vxi = 0. By applying these equations to the hehcity 
sums, we find 

^Q = ^{mQ + Ifh,. (22) 

For the K*, we replace the 75 with ^, and let tuq ms- The last substitution is strictly 
not rigorous, because we are not using a peaking approximation for the K* distribution 
amplitude. This will induce errors of order mx*/MB in the final results. Required for the 
calculation is the conjugate of the state, ^k* = 7°^lf*7°, 

^K* = -^^r{mK^+If). (23) 

Our calculation will include the exchange of a single hard gluon between the initial heavy 
or s quark and its spectator quark. The two diagrams are shown in figure 2(a,b), representing 
respectively the single gluon exchange between the heavy or s quark and its spectator. Here 
F represents the operator responsible for the heavy to light transition, as in the previous 
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(l-^)p ^(1-2/)/ {l-x)p ^ -^{l-y)p' 

(a) (b) 

FIG. 2. The two Born amplitudes contributing to Mq — > K* decay. The symbol F denotes an 
arbitrary Dirac matrix associated with the decay. 

section. These diagrams yield the contribution to the hard scattering amplitude Tp from 
perturbative QCD. We use the identity T"T" = |Jc where Ic is the 3x3 identity matrix of 
color space. The Brodsky-Lepage formalism tells us that the full amplitude at large w is 
given by 

< K*\sTQ\Mq >= J dxdyWK*{y)Tr{y,x)^Qix) (24) 

Both diagrams are infared (IR) divergent when the momentum flowing through the gluon 
line vanishes. Fig. (2b) also has an IR divergence associated with the strange quark going 
on-shell for small but finite gluon momentum transfer. These IR divergences are fictitious 
|T7| and irrelevant at large recoil. The region of y integration associated with the IR problem 



corresponds to the Drell-Yan-West |2^ end point region, which is assumed to be cut off by 



a Sudakov form factor. Physically, the mutually canceling effects of multi-gluon exchange 
cut off this region, due to the fact that the color singlet meson decouples from soft gluons 
in the IR limit. 

Following the cited literature, we apply a transverse momentum cutoff to suppress IR 
effects. In the rest frame of the heavy quark, we have the following picture. The momentum 
of the K* is given by mx*(cosh6', 0, 0, sinh6'). In this frame, the variable w = cosh^, and 
the spacelike gluon momentum is fcs = (1 — y)mK* sinh^^. This is the transverse momentum 
of the gluon, associated with transverse momentum of p — p'. Large transverse momentum 
means that this is larger than the typical momentum associated with the light degrees of 
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freedom, A. Therefore we apply the cutoff 

^-y> — r^^- (25) 

mK*VW — 1 

The calculation is done in Feynman gauge. To extract the heavy - light form factors, 
an expansion of the perturbative expression in powers of l/mq must be done. The leading 
terms for the gluon, strange quark, and heavy quark denominators (fc^ — in each case), 
in the peaking approximation, are given by 

Dg= {1- yfm\, + - 2Am/^*(l - y)w , (26a) 
Ds = m\, - + A^ - 2AmK*w = m^^. , (26b) 
Dq = 2mQ{A - {1 - y)mK*w) = itiqAq (26c) 

By comparing the perturbative expression to the heavy - light parametrization (^), we find 



= k{Iq - ^^^Is) , (27a) 

rriK* 

62 = -nilQ + Is) , (27b) 
2A 

^3 = -K Is , (27c) 

TTlK* 

^4 = 0, (27d) 
to leading order in the l/mq expansion. The constant k, is 



From general heavy quark symmetry arguments and the perturbative calculation done 
here, one finds that fq ~ Xj ^mq for large mass mQ, and therefore k is finite and non- 
vanishing in the infinite heavy quark limit. The scale of is set by the Mq meson mass. 
We will use for the B meson ag = 0.20. The functions Iq and Ig denote the integrals over 
the momentum fraction y from figure (2a) and (2b) respectively. 



l-e 



Iq= r ^dy^Uy)Pjr , (29a) 

Jo ^QJ-'g 
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W 

FIG. 3. The two integrals Fg (black curve) and Fq (gray curve) as defined in the text. The 
solid (dashed) curves are for the </>i {(pczz) K* distribution amplitude. 

where e denotes the IR cutoff given by eqn. (|25|) . We provide the exact integrals for arbitrary 
wave function and cutoff in appendix A. For the distribution amphtudes 0i and (p^zz defined 
above, they are plotted in fig. 3. From these figures, one can observe that the Drell-Yan- 
West region is turning over Jq at w ~ 2; that is, the loss of integral due to the IR cutoff 
is "beating" the IR divergence around this region. The integral Iq turns over at a smaller 
value of w because of the additional strange quark pole. For the remainder of this paper, 
we will study only the 0i case. 

To study the region (in w) of validity of this calculation, one can compare the asymptotic 
expansion, in powers of to the full result for Oi and 6*2. To leading orders in this 
expansion, 

ei{w) ~ , (30a) 
\nw 0.23 

02{w) ~ — . (30b) 

We note that 6i also possesses a. term but in our case it is numerically small. In 

fig. 4, the full and asymptotic (to order 1/w'^) values of 6i/k and 62/ n are plotted, for the 
01 distribution amplitude. For w < 2, less than two thirds of the integration region over 
y is included by the Drell-Yan-West cutoff, indicating a large soft cutoff dependence to the 

16 



2 3 4 5 6 7 8 

w 

FIG. 4. The un- normalized heavy- light form factors 9i/k, (gray) and O2/K (black), with their 
asymptotic behavior (dashed). 

amplitude. Our perturbative QCD result does not give a reliable estimate of the form factor 
in this regime. For w = 3, the Drell-Yan-West region cuts off about 18% of the momentum 
fraction integral, and the match between the full and asymptotic functions is reasonable, 
as shown in fig. 4. For w > Q the next to leading order w mK*/iTiQ corrections begin to 
dominate, and the leading order 6i do not meaningfully describe the amplitude. We will 
discuss the effect of the next to leading order corrections in sec. 5. 

Note that the parametric form of the thetas given by eqn. (pTl) has the correct heavy 
- heavy limit. As is taken to infinity, the integrals Iq and Ig remain finite and non- 
vanishing, and ^ = ^1 — is the leading contribution. 

IV. INTERPOLATING FUNCTIONS AND DETERMINATION OF 1^*^! 

For w less than or equal to about 2, interpolating functions for the heavy - light form 
factors are definitely required. Our strategy is to obtain such an interpolating function that 
is consistent with the data for D K*ev aX w = 1.29, matches the QCD calculation at 
w = 3.05, and is consistent with vector dominance ideas. For the calculation of Vts, we need 
to find interpolating functions for 61 and 62- As input to determine their forms, we use the 
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two data points at w = Wd, and the three ratios 

These three ratios are determined by perturbative QCD in the context of the Brodsky - 
Lapage formahsm. For the 0i distribution amphtude, ri = 0.2776, r2 = 0.4403, and = 
—0.6991. These values can be systematically improved by higher order (in a^) calculations 
|. As discussed in the next section, we have reason to believe that the contribution 
is a good approximation to the full result. We will consider a five parameter fit to the two 
form factors. Vector dominance ideas, that the amplitude is dominated by intermediate bs 
resonances, dictate that reasonable forms for the interpolating form factors are 1/w (single 
pole), l/w"^ (dipole), or a constant (subtraction). A subtraction is not considered for 02 
because the physical axial- vector form factor E {6i — W62), so that a constant term in 62 
would correspond to E w + . . . which is "inconsistent" with vector dominance. 
The interpolating functions used are 

e,=a(^ + -]+b, (32a) 

w J 

02 = a(^ + ^] . (32b) 
w J 

The theoretical constraints, that the slopes of 61, 62 and the ration 61/62 at wb match the 
perturbative result give 

c = -(ra - b/a)w% + nr^w^ , (33a) 

d = 2wB{r-i — b/a) — rir^w^ , (33b) 

e = + r2W% , (33c) 

/ = 2wb - r2wl 5 (33d) 

so that experiment fixes a and b. The fit to experiment is shown in fig. 5. The dashed 
lines refer to the parametrized fit, and solid lines refer to the full IR sensitive functions. 
By construction, the fit and original function match very well at w = wb- (Actually, the 
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FIG. 5. The normalized heavy- light form factors 6i (gray) and 02 (black), with their 
parametrizations with respect to the K* data (dashed). 

full perturbative functions extrapolate to the data reasonably well without any fits. This 
is another indication that the perturbative calculation is sensible.). The propagation of 
experimental uncertainty for V and Ai from wd given by eqn. (5) to wb via our interpolating 
functions is shown in fig. 6. The dashed lines denote the boundaries of one standard deviation 
envelopes about the mean values. 

The relevant function for the calculation of |Vts| is the combination of tensor form factors 
1^ as discussed in sec. 1, and defined in sec. 2. From eqns. (|13D and the interpolating 
function displayed in fig. 6, we find 

I^W + 5(^)|(g2 = 0) = 0.78 ± 0.16 (21%) . (34) 

By combining the CLEO data for B K*'~f, the B lifetime, and the theoretical decay rate, 
given in sec. 1, with this result, we find 

\Vts\ = 0.035 ± 0.010 (28%) . (35) 

Here we have quoted only the experimental uncertainty, about 21% from D — > K*ii' data, 
and 19% from B — > K*'~f data, which in turn are added in quadrature. The theoretical 
uncertainties are discussed in the next section. 
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FIG. 6. The parametrized heavy-light form factors 9i (gray) and 62 (black), with one standard 
deviation of uncertainty (with respect to the D K*iv data 1^(0) and A1(0)) denoted by the 
dashed curves. 

This result compares quite well with the standard model value \Vts\ = .041 ± 
.007 (17% total uncertainty) extracted from Vcb and the unitarity of the CKM matrix (dis- 
cussed in sec. 1). Therefore our analysis shows that the B —* K*'~f data agrees with the 
standard model, and should be used to place constraints on non-standard model physics. 



V. THEORETICAL UNCERTAINTIES AND Fb 



Now consider theoretical corrections to our result for \ Vts \ given by eqn. ([351) . We consider 
corrections to the heavy quark symmetry relations, to the as{mB) matching conditions, and 
uncertainty in the top quark mass. 

The corrections to heavy quark symmetry relations are due to physics which violates 
the assumption that in the rest frame of the heavy meson, the heavy quark is also at rest. 
The first of this type is due to the QCD interactions with gluons and light quarks in the 
heavy meson. The correction factor is [p!0| , p6| A/2mg, where A is the mass difference between 
the heavy meson and heavy quark. For a conservatively low value of rric = 1.3 GeV, this 
correction is roughly 15% for the D meson. 
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FIG. 7. The percentage of heavy quark symmetry violating w mK*/mB corrections for 9i 
(dashed gray), 62 (dashed black), and A + B for the B K*j process (solid black). 

The second type of HQET correction occurs when the momentum transfer between initial 
and final mesons is large. Then, the emission of a hard gluon by the heavy quark before 
its decay can give the heavy quark a large velocity with respect to the meson rest frame. 
The naive estimate of this effect is [|TU|,|TT],^ {v ■ p')/2mq ^ 1/4 at zero q^. However, 
studies of these QCD corrections indicate that the problem caused by the leading gluon 
exchange diagram is strongly suppressed by an order of magnitude 0,^. We have used the 
perturbative QCD calculation of sec. 3 to estimate these corrections. They are of the form 



w rriK* 
rriB 



C 



(36) 



where 



C 



1- , <I^Hy) 

'^A^o 



(37) 



This correction comes from from the diagram of fig. (2a). We evaluated it for = 0i as 
defined in sec. 3, and plot the the ratios AOi/Oi, A6'2/6'2, and A{A + B)/{A + B) in fig. 7. 
As w grows much larger than m^^. /jt^b, these corrections actually dominate the asymptotic 
behavior. The solid line of fig. 7 evaluated aX w = wb gives an estimate of these corrections 
to \Vts\ of about 12%. While the wd mK*/rnD corrections are impossible to estimate via 
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perturbation theory because wd so close to one, heavy quark symmetry dictates that they 
should have the same form as the wb itlk'/itlb corrections, that is, they are systematically 
correlated. The experimental values for the ^'s at wd includes these corrections, while our 
perturbative analysis does not. Therefore these errors systematically compensate each other 
in the final result, so we use 12% uncertainty total, including both corrections. 

Another source of HQET violations are perturbative correction which cannot be included 
in the definitions of the heavy - light form factors 6i. From our perturbative calculation, 
we find that the dominant source of these corrections, which are proportional to (1 — ;^)/2, 
are precisely the same algebraic form as the w ■ mK*/rnB corrections to the 6i as discussed 
above. It is clearly difficult to assign a percentage of uncertainty to these corrections, since 
we have not calculated the precise way in which they feed into our calculation of \Vts\, so we 
make a naive estimate of 6% total uncertainty from them. This value is smaller than the 
naive estimate of wrriK* /ttib ~ 1/4, but not by an order of magnitude as was hoped P,p7|. 

We also need to estimate corrections to our perturbative QCD calculation of the matching 
conditions given by eqns. (^). That is, we wish to estimate the uncertainty due to the 
running of the form factors from w = wj:) to w = wb- It has been suggested in the literature 
that the ag contribution to the form factors may be producing only 10% of the full result 
0]. In this case, the matching conditions would be meaningless. 



To check this, we use eqn. (^) to determine /b- From the D K*iiy data and our fits, 
we find 



^ = fK'fB ^ 3 = O-O'^S ± .022(28%) , (38) 

mi-* 9 



^JmK*fnB Anas 

where we have quoted experimental uncertainties. To extract from eqn. 
«s(^b) = 0.20, and fx* = 212 MeV from QCD sum rules p4[. This sum rule estimate of 
fx* is supposedly is good to about 20%. This yields 

fB ~ 420 MeV. (39) 

This is a rather large value of /b, when compared to numerous lattice estimates of 200 to 
300 MeV. However, it is not an order of magnitude larger than these values, which would be 
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the case if the corrections contributed only 10% to the full result for the form factors. In 
addition, the calculation of fs done here assumes the leading HQET relations between the 
D and B meson systems. Decay constants are well known to receive large 1 / itlq corrections 
in the HQET |]10[, and these should be taken into account before a real prediction of fs 
is made via D K*lv and B K*'y data. A final note on this topic is that there is 
new evidence that the lattice calculations may systematically underestimate heavy quark 
decay constants. The recently measured fo^ decay constant's central value is about 1.5 



times the central values of the lattice estimates. We conclude that the order as calculation 
is indeed dominating the form factor calculation at w = wb- And a more proper way of 
calculating fs would be to include all l/rriB corrections in the HQET. 

Since the normalization of the form factors is set by experiment in our case, the error 
in Jb is not directly correlated with error in \Vts\- Our input for \Vts\ from perturbative 



QCD are the ratios ri,r2,r3 defined in equation (0). We currently have no systematic way 
of estimating their uncertainties due to higher order corrections, other that to apply that 
standard Os/tt rule for the next to leading order correction. This yields a naively small 
error estimate of 7%. Clearly, an explicit calculation of the type given by Field et. al. 
would pin this uncertainty down further. The uncertainty in as{mB) is about 10%. 



Our uncertainty in the soft physics which goes into the perturbative calculation can also be 
estimated. The peaking approximation for the heavy meson distribution amplitude yields 
about A/rriB = 10%, and the uncertainty in the K* wavefunction about mx*/TnB = 17%. 

In addition, our perturbative QCD calculation is sensitive to the IR cutoff used to regulate 
our perturbative momentum fraction integrals. This is the standard "solution" to the IR 
problem in heavy - light systems as discussed in the literature 0,^,0. This aspect of the 



calculation definitively needs to be improved to given more confidence to the application 
of perturbative QCD to heavy quark systems. (Our philosophy in this paper has been 
to apply this "well known" QCD technology.) One way of estimating uncertainty due to 
this ambiguity is to perform the calculation for various distribution amplitudes. While the 
results were not displayed, we found that the (pczz distribution amplitude gives essentially 
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the same result for \Vts\, to within about 10% uncertainty. The similarity of the result is 
already evident from fig. 3. Finally, theoretical uncertainty from the top quark mass adds 
about 8%. 

We combine the theoretical uncertainty into two parts. The first part consists of uncer- 
tainty due to corrections in HQET and the top quark mass. They give about 22% when 
added in quadrature, with the dominant contributions comming from w ■ rnK'/friQ type cor- 
rections. The second type of theoretical uncertainties come from running the form factors 
from wd to wb via the perturbative QCD matching conditions. We estimate uncertainty 
of about 23% from perturbative QCD uncertainties. Clearly, these are very difficult to es- 
timate because we don't know how the QCD corrections feed into the matching conditions 
eqns. (|3TD. However, because we are using perturbative QCD rather than simply making 
a pole or dipole ansatz for the form factors, we can make an estimate that can be system- 
atically improved as the perturbative calculations become more sophisticated. Hence we 
find a total estimated theoretical uncertainty of about 32% from all sources when added in 
quadrature. 



VI. SUMMARY AND CONCLUSIONS 

We find a value for the mixing angle of about 

\Vts\ = 0.035 ± .010 ± .011(28% + 32%) (40) 

where the second one a uncertainty is a naive estimate of total theoretical errors. It is naive 
because it is an estimate of higher order corrections that we have not calculated. We believe 
that this factor can be substantially reduced by further work on l/mg HQET corrections 
to the heavy quark matching between D and B systems, and an order calculation of the 
perturbative QCD matching conditions for the heavy - light form factors a.t w = wb- The 
most disturbing theoretical uncertainty seems, at this point, to be the IR sensitivity of the 
perturbative QCD calculation, which we have estimated to be 10%. 
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Our value compares very well with the standard model result of \Vts\ = 0.041 ± 
.007(17%), that is, our result is consistent with three generations of quarks, and the standard 
model contributions to the 6 — > 57 penguin. 

We can compare our result with recent lattice estimates of the hadronic tensor form 
factor for B K*-f decay. The UKQCD group estimates 2Ti = \A + B\ = .58ttl, and 
with an additional assumption of spectator mass independence, |A + i?| = .60^^°. Bernard 
et. al. H estimate \A + B\ = AOt°u 

using a pole form ansatz. These results should be 
compared to our result of |yl + i?| = 0.78 ± .16 ± .25 (21 ± 32)%, where we quote separately 
experimental and theoretical uncertainty. (There is less experimental uncertainty than for 
\Vts\ since we do not require B K*'~f data for this parameter). 

The QCD sum rule technique has also been applied to help determine the B K*'y form 
factors ||3^, and a single pole ansatz, with no subtractions, has been recently used |0. As 
discussed in sec. 2, this ansatz is hard to reconcile with leading order HQET, so that A/rriB 
corrections play a significant role in determining the form factors by this method. However, 
both of these results are in rough agreement with ours. 

As part of a self consistency check on the perturbative part of our analysis, we find the B 
meson decay constant fs ~ 420 MeV with 28% ±32% ±20% uncertainty, where the last 20% 
is the generic uncertainty in the fx* decay constant from QCD sum rules. While our value 
of \Vts\ is not directly correlated to this value, since we fix the normalization of our form 
factors from experiment, this value does indicate that the perturbative QCD calculation is 
correctly estimating the exclusive decay rate at large meson - meson recoil. 
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APPENDIX A: CALCULATION OF INTEGRALS Iq AND Is 

The QCD calculation for 9i, 62 and ^3 is written in terms of two integrals over the light 
quark momentum fraction y. These integrals Iq and Ig arise from the Feynman diagrams 
containing a heavy or strange quark propagator, figure 2a and 2b respectively. They are 
given by: 



Ic 



^ dy(p^K*{y)- 



m 



K* 



^QDg 

'■1 m?. 



where Aq, Dg and Dg are given in eqns. (26) and y — 1 — y. It is convenient to rewrite these 
integrals in terms of 

h = -2wIq = C dy 1 ^—-2 1 (41) 



and 



rriK* " Je " y'^ + 2^^wy 

By partial fractions these integrals can be written in terms of a single function 

h = ^(/(«i) - /(«3)) , (43a) 

ai — as ai — 02 oi — 02 

l2^^—ifia2)-fias)) , (43b) 
a2 — as 



where 



ai = , (44a) 

mK*w 

A 



02 = (w + - 1) , (44b) 

as = -^(w - Vw"^ - 1) , (44c) 
rriK* 

26 



gives the position of possible poles in 1 — y, (before the Drell-Yan-West region has been 
cutoff), and the function /(a^) is defined by 



f{a^) = I ' dy$^ . (45) 



The value of /(oj) involves the cutoff prescription e and the form of the K* wave function. 
In this paper we have employed the w dependent cutoff e — ^ ^* 
function is taken to be ^^^.(y) = 6y(l — y), then 

/(cj) = 3 — 6ai + 3e^ + 6{ai — l)e + 6ai(l — a^) In |1 — ai\ — 6ai(l — a^) In |e — ajj (46) 

Note that this function, together with 7q, Ig and ^i, ^2, ^3, depends only on and w, and 
is independent of uiq as required by heavy quark symmetry The explicit form for the 9i 
involves logarithms the largest of which goes like ^ as found in the asymptotic expansions 
given by eqn. (30). 
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